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Based on a series of lattice calculations we determine the ratio Fx / -FV 
in QCD. With experimental data from kaon decay and nuclear double 
beta decay, we obtain a precise determination of \Vus\. Our simulation 
includes 2 -f 1 flavours of sea quarks, with three lattice spacings, large 
volumes and a simulated pion mass reaching down to about 190 MeV for 
a full control over the systematic uncertainties. 
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1 Introduction 



The measurement of CKM matrix elements involves examining weak decays of hadrons. 
The theoretical description of these decays involves computing strongly interacting 
effects of quarks and gluons inside the hadrons. These effects should be determined 
in a fully non-perturbative way. Any weak decay can be written as the product of 
three factors 

Weak decay = {Kinematic} x {CKM element} x {Non-perturbative QCD}. (1) 

Computing the non-perturbative QCD effects is a difficult task. A first princi- 
ple analytical approach does not yet exists. We have to live either with numerical 
simulations or with models of QCD. 

Here we will use lattice QCD as a numerical tool to compute non perturbative 
QCD meson decay contants. 

1.1 Light meson decay constants and \Vus\ 

We will follow a proposal by Marciano [1] of examining the ratio of decay widths 
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where this ratio, even including electromagnetic corrections, are experimentally know 
with an 0.4% precision. Masses are known at least with a relative precision of 3 x 
10~^. This implies that an accurate determination of the ratio FkIF.,^ can be used to 
determine |V^s|/|V^d|. The most recent update of the Flavianet kaon working group 
is 12]: 

|Ks| Fk 
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0.27599(59) . (3) 
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This may be used with = 0.97425(22) [3] from super-allowed nuclear /5 decays 
and the lattice determination of Fk/Ft^ to obtain \Vus\- 



2 Simulation details 

We simulate 2 degenerate light clover fermions representing the u and d quarks in the 
strong isospin limit plus a third heavier clover fermion that represents the strange 
quark. Whereas the value of the strange quark mass is held fixed close to its phys- 
ical value light quark masses are varied through a significant range so that our vr 
mesons have masses in the range 190 MeV<M^<460 MeV. This allows a controlled 
extrapolation to the physical point M^r ~ 135 MeV (see Figure [Ta|) . 

We simulate at three different lattice spacings a ~ 0.124 fm (/3 = 3.3), a ~ 0.083 fm 
(/3 = 3.57), and a ~ 0.065 fm (/3 = 3.7) for a controlled extrapolation to the continuum 
limit. 
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(a) Overview of our simulation points in the (b) One of the 1512 fits that leads to our final 
plane (M^, y/2M]^ - M^). determination of Fk / 

Figure 1: Overview of simulation points and a representative fit. 



All our simulations are done in large volumes M7r-L>4, so that finite volume cor- 
rections are below 1%. 

More details about both our fermion and gauge actions can be found in [1]. 



3 Systematic uncertainties. 

To obtain a QCD results from lattice data one has to perform several extrapolations: 
one has to extrapolate data to the physical value of the quark masses (chiral extrap- 
olation), one has to to the infinite volume limit, and to the continuum limit. One 
has also to deal with the contribution from exited states and the possible uncertainty 
associated with determining the QCD coupling (i.e. setting the scale). All these 
considerations should be taken into account and each of them should be reflected in 
the final uncertainty of our result. 

A detailed description of the method sketched here can be found in the original 
work [5]. 

3.1 Scale setting issues 

2-1-1 flavour QCD has three independent parameters: the light quark masses, the 
strange quark mass, and the gauge coupling. One needs three experimental inputs to 
fix these parameters in our simulations. 

To fix quark masses one uses the experimental values of the masses of the vr and 
K mesons, M^^ and M^- 

To set the lattice spacing, one uses a physical state whose mass does not strongly 
depend on quark masses. We used both the Vl and H states, the differences between 
these two choices measures the error associated with our choice of physical input. 

All physical inputs are corrected to take into account that we simulate in the 
isospin symmetric limit without electromagnetic interactions. 
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3.2 Chiral extrapolation 



There exist basically two approaches that can guide us in the extrapolation to the 
physical point. We can use an effective field theory approach (xPT), either in the 
two or three flavour setups. 

Another possibility comes from the observation that physical quantities like decay 
constants or masses have an analytical expansion in the quark masses as long as this 
expansion performed around a regular point far from the singular point rriq = 0. If our 
lattice data is close enough to the physical point these analytical expansions should 
be a fair approximation to the behaviour of the physical quantity up to the physical 
point. 

We will flt our data to a total of 7 different expressions that correspond to one of 
the previously mentioned approaches. 

For the case of three flavour xPT our data is well fltted using the NLO relation for 
the ratio Fk/Fj^ [6]. In fact up to higher order terms, the NLO relation can be written 
in three different ways. The difference between these three expressions measures the 
contribution of higher order contributions. 

The ratio Fk/F^^ can be obtained from the pion mass dependence of Fk in the 
heavy kaon two flavour xPT formalism and the pion mass dependence of F^^ in SU{2) 
xPT [7J. Our data is well fitted by the NLO relation, and up to higher order terms, 
we have two different ways of expanding the ratio of decay constants. The difference 
between these two expressions are higher order contributions. 

Finally the ratio of decay constants have a perfectly regular expansion, as soon as 
this expansion is done far from the singular point ruq = 0. In fact our lattice data is 
well fitted by a simple polynomial expansion around a regular point that lies midway 
between the physical point and our heaviest pion ensemble. One can also fit the data 
to a rational approximation, that can be seen as a Fade like ansatz. This gives a total 
of 2 expressions based on an analytical expansion around a regular point that, again, 
differ in higher order terms. 

It is important to remark that the total of 7 formulae differ both in the approach 
(three, two flavour xFT or analytical expansions), and in the higher order contri- 
butions within a framework (NNLO contributions in xPT or higher powers in the 
analytical expansions.) 

3.3 Continuum limit 

The ratio Fk/Ft^ is a flavour breaking term. This means that no matter the value of 
the lattice spacing we will always have Fk/F^^ = 1 when niud = ms- This suggests 
that cutoff effects should be oc (M^ - M|). 

Our action is formally only improved up to order 0{asa), but we have numerical 
evidence |1] that cutoff effects scale as C(a^). 

In particular for the present case, cutoff effects are small. This can be understood 
by looking at SU{3) xPT, where the leading correction cancels for the ratio of decay 
constants. For our particular dataset they are statistically consistent with zero (i.e. 
the statistical uncertainty of our data are bigger than lattice spacing corrections.) 
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Taking these considerations into account we choose to parametrise cutoff effects 
in our data by adding or not a term 

r a(Mi - M2)//iQCD 
= c X <^ or (4) 
[ a\Ml - Ml) 

to any of the previous 7 formulae that gives the dependence of Fi^jF^^ with meson 
masses. 

3.4 Finite volume corrections 

For large enough lattices, finite volumes effects scale as oc e~^^^^. All our simulation 
points respect the bound M^>4, making finite volume corrections small. In fact our 
data is well fitted without any kind of finite volume correction, indicating that the 
size of the correction is smaller than the statistical accuracy of our data. 

Nevertheless, corrections due to finite volume are known in xPT. Meson masses 
and decay constants receive corrections that have been computed. We decided to 
correct the decay constants with a 2 loop expression [8]. To estimate the uncertainty 
of this correction, we use the 1 loop expression together with an upper bound on the 
correction (coming from correcting only F^). 

3.5 Excited state contributions 

Masses and decay constants are obtained from the large time behaviour of correlators, 
i.e. 

m - {^'i:m.u](x)[n-„m) °g^^ <°fe"i'^"(g)><^"(g)i"^'''i°) e-''-- 

(5) 

We choose for each value of the couphng /3 the start of the fitting interval t^iin 
so that it is strongly dominated by the ground state. There are several reasonable 
possibilities for this, and we choose to use 18 possibilities that will receive different 
tiny contributions from excited states. These different fitting intervals will be used 
to estimate the uncertainties due to excited state. 

4 Results 

The procedure described above lads to a total of 2 x 18 x 7 x 2 x 3 = 1512 different 
ways of obtaining the ratio Fk / F^^ in the continuum, at the physical mass point and 
in the infinite volume limit (see Figure ITbl) 

Since not all the previous procedures describes our data equally well (although 
the x^/dof for all the fits is close to 1), we weight each of the values with the quality 
of the fit to produce a distribution of values. 
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The median, typical result of our analysis, is taken as our final result. The 16- 
th/84-th percentiles, that measures the "spread" of values, give our final systematic 
uncertainty. To determine the statistical accuracy of our results, we resample our 
configurations 2000 times, and the standard deviation of the median over these 2000 
samples is taken as the statistical uncertainty. 

Our final result is 

^ = l-192(7)3tat(6)syst (6) 
i'n phys 

Using the value |Kd| = 0.97425(22) [3J we obtain 

iKsl = 0.2256(18). (7) 

Using also \Vub\ = (3.93 ± 0.36) x 10"'^ the first row CKM unitarity relation is 

iKdP + iKsP + \Vub\^ = 1.0001(9) (8) 

with no signal for physics beyond the standard model. For details on the calculations 
see [S]. 
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